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When  low-energy  ions  are  incident  on  a  solid  surface,  a  large  frac¬ 
tion  is  neutralized  via  an  Auger  or  resonance  process.  Some  useful  in¬ 
formation  on  details  of  the  neutralization  process  can  be  obtained  by 
studying  the  ion  fraction  of  backscattered  particles.  An  accurate  evalu 
ation  of  the  characteristic  velocity,  which  determines  how  fast  the  ion 
fraction  varies  with  the  incident  ion  velocity,  is  expected  to  be  partic 
ularly  useful.  In  this  paper  we  give  a  theoretical  analysis  of  experi¬ 
mental  ion  yield  data  that  exist  in  the  literature.  It  is  shown  that  a 
successful  interpretation  of  the  ion  yield  data  can  be  made  by  taking 
into  account  the  variation  of  the  characteristic  velocity  with  respect 
to  the  incident  ion  energy. 


1 


Introduction 


Much  experimental  and  theoretical  work  on  neutralization,  excitation 
and  deexcitation  of  ions  at  a  solid  surface  has  been  reported  in  the  lit¬ 
erature  11-7].  Theoretically,  low-energy  (£  5  keV)  surface  processes 
have  successfully  been  described  by  the  exponential  model  of  Hagstrum 
[8-10]  in  which  the  rate  of  surface  processes  is  assumed  to  decrease  ex¬ 
ponentially  with  the  ion-surface  separation.  Experimental  electron  and 
photon  yield  data,  in  particular,  have  given  strong  support  to  the  model. 
Recently,  however,  a  nonlinear  dependence  of  the  logarithm  of  the  ion 
yield  upon  the  inverse  of  the  initial  (or  final)  ion  velocity  has  been 
observed  by  several  groups  [11-14],  in  contradiction  to  the  prediction 
of  the  model.  This  is  generally  believed  to  be  an  indication  that  ions 
are  subject  to  more  than  one  neutralization  mechanism.  Theoretical  stud¬ 
ies  in  the  past  have  shown  that  reasonable  agreement  can  still  be  ob¬ 
tained  if  one  considers,  in  addition  to  surface  neutralization  (Auger  or 
resonance  neutralization),  neutralization  and  ionization  resulting  from 
atomic  collision  between  the  incident  ion  and  a  surface  atom.  Thus,  the 
ion  trajectory  is  conveniently  divided  into  three  parts  [15-17]:  the  in¬ 
coming  trajectory,  the  region  of  violent  atomic  collision  and  the  out¬ 
going  trajectory.  Along  the  incoming  and  outgoing  trajectories,  Auger 
or  resonance  neutralization  may  occur  while  collisional  neutralization 
and  ionization  of  the  atomic  type  take  place  in  the  region  of  violent 
collision. 

In  our  earlier  work  [18]  we  have  introduced  a  model  which  presented 
a  different  view  of  the  ion-surface  interaction.  According  to  the  model, 
the  observed  nonlinear  behavior  of  the  ion  yield  does  not  require  the  in¬ 
troduction  of  atomic  neutralization  and  ionization.  The  idea  stems  from 


Che  observation  that  the  "characteristic  velocity"  (which  determines  how 
fast  the  scattered  ion  yield  decreases  with  the  ion  velocity)  is  itself 
a  function  of  the  ion  energy,  resulting  in  a  more  complex  dependence  of 
the  logarithm  of  the  ion  yield  on  the  inverse  of  the  ion  velocity  than  a 
simple  linear  relationship.  Good  to  excellent  agreement  has  been  ob¬ 
tained  in  the  preliminary  comparison  of  our  model  with  experimental  ion 
and  photon  yield  data. 

In  this  paper  we  further  develop  the  model  introduced  in  ref.  18  so 
that  it  can  be  applied  to  a  wider  class  of  experimental  data  on  nonspec- 
ularly  as  well  as  specularly  reflected  particles.  We  then  report  the 
result  of  an  extensive  test  of  our  model  against  experimental  ion  yield 
data  on  various  ion-surface  combinations.  Since  for  the  present  work  we 
are  interested  only  in  the  charge  distribution  of  backscattered  parti¬ 
cles,  there  is  no  need  Co  consider  excitation  and  deexcitation  processes 
which  contribute  negligibly  to  the  ion  yield  at  low  energies  of  our  inter¬ 
est.  Our  analysis  is  centered  on  an  accurate  determination  of  the  char¬ 
acteristic  velocity,  necessary  for  the  calculation  of  the  ion  yield.  It 
is  our  hope  that  the  present  work  gives  useful  insight  into  the  complex 
problem  of  ion-surface  interaction  at  the  fundamental  level.  From  the 
practical  viewpoint,  the  accurate  determination  of  the  characteristic  ve¬ 
locity  should  be  useful  for  surface  analysis  because  the  characteristic 
velocity  is  sensitive  to  details  of  the  surface  such  as  the  molecular 
composition,  electronic  state  and  surface  contamination. 

In  section  2  we  give  a  general  description  of  our  model.  In  partic¬ 


ular,  we  present  a  detailed  analysis  of  the  characteristic  velocity  for 
both  specularly  and  nonspecularly  reflected  particles.  Extensive  quanti¬ 
tative  comparison  of  the  prediction  of  our  model  with  experimental  data 


that  exist  La  the  literature  Is  presented  in  section  3.  Finally 
cussion  is  given  in  section  4. 


A  detailed  description  of  our  model  in  terms  of  rate  equations  has 
already  been  given  in  ref.  18.  Since  for  the  present  discussion  we  are 
not  interested  in  photon  yields  or  other  data  related  to  excitation  and 
deexcitation  of  ions,  surface  neutralization  (Auger  or  resonance  neutral 
ization)  is  the  only  process  we  need  to  consider.  We  also  note  that 
atomic  neutralization  and  ionization  are  not  allowed  in  our  model.  The 
rate  equation  then  becomes 

dNL(t) 

— -  -r(t)NI(t),  (1) 

dN2(t) 

dt  •  TCON^t),  (2) 

where  is  the  number  of  ions,  the  number  of  neutral  atoms,  and  r  is 
the  neutralization  rate  assumed  to  decrease  exponentially  with  respect 
to  the  ion-surface  separation  z, 

T(z)  =  A  exp(-az).  (3) 


Assuming  a  straight-line,  constant-velocity  trajectory,  i.e.,  assuming 


z-z 


t  <  0  , 
t  ,  t  >  0  , 


-vuc  * 


(4) 


where  z  is  the  distance  of  closest  approach,  and  v.  and  v,  are,  re- 

O  la.  fa.  * 

spectively,  the  normal  components  of  the  initial  (v^)  and  final  (v^) 
velocity  of  the  ion,  T  can  be.  expressed  in  terras  of  time  t  as 


av  exp(av,  t)  ,  t  <  0  , 
C  l_i_ 


avc  exp  (-av^t)  ,  c  >  0 


Here  v  is  the  characteristic  velocity  given  by 
c 


vc  =  -  exp  <-azo)  , 


and,  assuming  a  single  binary  collision,  v^  and  v^  are  related  by 


v £  **  u(9)vi  , 


where 


u(0)  ■  — tbj  cos0  +  -  sin20j  , 


0  is  the  scattering  angle,  and  ra^  and  m2  are  the  masses  of  the  incident 
ion  and  the  surface  atom,  respectively* 

We  consider  the  case  when  Che  incident  particles  are  singly  charged 
ions,  corresponding  to  the  initial  condition,  ■  N,  *  0.  Eqs 

(1)  and  (2)  with  this  initial  condition  can  immediately  be  solved  to  yield 


—  •« 

N,(»)  =*  N  exp  -v  (——  +  , 

1  K  vtJ 

No(0#)  =  N  -  N  exp  -v  ^ 

2V  H  c\vix  vfj_ ) 


The  ion  fraction  f  is  then  given  by 


(9) 


Ni(») 

Nj  («)+N2(°°) 


cWj.  vf  j.  / 


In  the  standard  exponential  model  [8-10],  the  characteristic  velocity  vc 
is  assumed  independent  of  ion  energy.  Thus,  when  log  f  is  plotted 
against  +  ^  >  one  would  obtain  a  straight  line  and  the  character¬ 

istic  velocity  could  be  determined  by  the  slope  of  the  line.  We  assert, 
however,  that  v^  is  not  independent  of  energy  as  can  clearly  be  seen 
from  eq.  (6):  the  characteristic  velocity  depends  on  Zq  which  in  turn  var¬ 
ies  with  the  ion  energy.  Hence,  according  to  our  model,  log  f+  will  not 

vary  linearly  with  ( — —  +  — — ^  .  An  accurate  determination  of  the  charac- 

W.L  "fa./ 

teristic  velocity  v^  and  its  comparison  with  experimental  data  thus  pro¬ 
vide  a  severe  test  of  our  model.  Below  we  give  a  detailed  analysis  of 
the  characteristic  velocity  and  discuss  its  dependence  on  the  ion  veloc¬ 
ity  and  other  collision  parameters. 

Let  us  first  consider  specularly  reflected  particles  »  0  -  ij»;  see 

fig.  la).  At  the  distance  of  closest  approach,  the  total  energy  avail- 

1  2  2 

able  to  the  ion  is  approximately  m^(v^  +  ).  Assuming  that  the  compo¬ 

nent  of  the  ion  velocity  parallel  to  the  surface  does  not  change  much 
during  collision,  we  have 


1  /  2  .  2.  1  ,2.2.  2  _  "bz0 

j  m^v  +  v{)  -  j  m^v  +  vf )  cos  <j>  +  Be 


In  eq.  (10)  we  have  also  assumed  that  the  ion-surface  interaction  at  the 
distance  of  closest  approach  is  the  Born  Mayer  type.  Solving  eq.  (10)  for 
zq,  we  obtain 


"22  2 
1  mi(vj+vf)  sin  ifi 

*  b  108  - 4B - 


Substituting  eq.  (11)  into  eq.  (6),  we  obtain  the  desired  expression  for 
v£  for  specularly  reflected  particles: 


v 

c 


A 

’mi(v^+vf)  sin^"|a^bM  A 

2  o  ^ 

■m]Vi(l+ufc)  sinnjr 

a 

4B  J  3 

4B 

a/b 


(12) 


If  we  neglect  the  elastic  energy  loss  of  the  ion,  i.c.,  if  v^  *  v^,  then 
eq.  (12)  reduces  to  the  formula  derived  in  ref.  18.  According  to  eq. 
(12),  the  characteristic  velocity  V£  is  an  increasing  function  of  the 
initial  (or  final)  ion  velocity  at  a  given  angle  of  incidence  i|»,  consist¬ 
ent  with  the  observation  of  Verhey  et  al  [17],  Bertrand  et  al  [12]  and 
MacDonald  et  al  [19,20].  Eq.  (12)  also  indicates  that  the  characteristic 
velocity  is  an  increasing  function  of  the  angle  of  incidence  at  a  given 
ion  energy,  consistent  with  the  observation  of  Bertrand  et  al  [12]  and 
Overbury  et  al  [14]. 

We  now  consider  nonspecular  reflection.  The  ion  which  is  nonspecu- 
larly  reflected  experiences  the  repulsive  force  mainly  along  the  symmetry 
axis  of  the  trajectory  (z'  axis  in  figs,  lb  and  lc).  The  ion-surface 
potential  in  this  case  is  therefore  assumed  to  be  given  by  V  *  Be  bz  in 
the  neighborhood  of  the  distance  of  closest  approach.  For  the  case 
i|>  <  0  -  <|>  (fig.  lb),  the  distance  of  closest  approach  is  arrived  at  be¬ 
fore  the  ion  reaches  the  middle  point  of  the  trajectory  at  which  most  of 
the  elastic  energy  transfer  is  assumed  to  occur.  We  therefore  assume 
that,  at  zq,  both  the  total  energy  available  to  the  ion  and  the  paralell 
component  of  the  ion  velocity  have  not  changed  much  from  their  initial 
values.  Thus,  we  may  write 


8 


'!«'3 

;*»  Jl 


Ml 


;*$l 


y  ralvi  *  ~2  mlVi  cos  'J'  +  "e  • 


Noting  that  the  z  axis  and  z'  axis  make  an  angle  of  ^  y  -  ip^  ,  eq.  (13) 
immediately  yields 


'8  A  ,22 
,2  “  V  ,  [ mlvi  sln  ♦ 

k -  lo£  (  - r5 - 


z  =  - 
o 


and  therefore 


v  =  — 
c  a 


/  2  •  2,\f  cos  -ip) 

A  /mlvi  sin  'l>y>  2  r 


The  argument  goes  similarly  for  the  case  ip  >  0  -  <p  (fig.  lc)  except 
that  the  ion  passes  the  middle  point  of  the  trajectory  before  it  comes 
closest  to  the  surface,  and  therefore  at  zq  the  total  energy  of  the  Lon 
and  the  parallel  component  of  the  ion  velocity  are  not  much  different  from 
their  final  values.  We  then  have 


1  21  2  2/q  ,  .  A  _  -bz' 

y  TIljVf  =  y  m^  cos  (0  -  lp)  +  Be  ° 


'  0  \ 
,*-y 


2  2 

jVfsin  (0-'p) 


m^Vf  sin^(0-ip)\  b  C08^  2^ 


■•w- v v ■.« .--v- -.v>vc 


Eqs.  (15)  and  (18)  show  that,  at  a  given  energy  and  scattering  angle,  the 
characteristic  velocity  v£  decreases  as  the  trajectory  moves  away  from 
the  specular  configuration,  consistent  with  the  observation  of  Bertrand 
et  al  [12].  One  can  also  observe  that,  at  a  given  ion  energy  and  inci¬ 
dent  angle  (or  exit  angle  if  <0  -  t]»),  slowly  increases  as  the  scat¬ 
tering  angle  9  is  increased.  In  the  next  section  a  detailed  quantitative 
comparison  of  the  analysis  given  here  with  experimental  data  is  presented 
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3.  Comparison  with  experimental  data 
3.1.  H+  -  C 

Hydrogen  ions  and  atoms  backscattered  from  a  solid  surface  generally 
show  a  broad  energy  distribution  peaked  at  a  velocity  lower  than  that 
given  by  eq.  (7a).  In  such  a  case  our  analysis  based  on  the  assumption 
of  a  single  binary  collision  should  apply  only  to  a  small  fraction  of 
particles  that  are  "kinematically"  Scattered,  i.e.,  to  those  that  are 
scattered  with  velocity  given  by  eq.  (7a).  Earlier,  Overbury  et  al  [14] 
considered  only  such  particles  in  their  study  of  H+-graphite  collision. 
From  their  experimental  ion  yield  data  they  calculated  the  ion  fraction 
of  kinematically  backscattered  particles  using  scattering  cross  sections 
calculated  for  a  screened  Coulomb  potential.  Our  analysis  therefore  can 
be  directly  compared  with  their  results. 

Our  analysis  of  the  1  keV  specular  data  of  Overbury  et  al  was  given 
in  ref.  18.  With  the  elastic  energy  transfer  taken  into  account,  we  ob- 

cl  A 

tained  slightly  larger  values  for  the  parameters  and  —  than  given 

in  ref.  18:  -  0.38  (compared  with  0.37  in  ref.  18),  —  a  1.2  x  10®  cm/sec 

b  a 

g 

(compared  with  1.1  x  10  cm/sec  in  ref.  18).  Here  we  first  show  the  re¬ 
sult  of  the  analysis  of  the  2  keV  specular  data  of  Overbury  et  al.  Shown 
in  table  1  are  their  calculated  values  of  the  ion  fraction  f+  and  the  cor¬ 
responding  values  of  v,  ,  taken  from  fig.  3  of  their  paper.  Also  shown 

t-L. 

are  the  values  of  <J>,  v.  and  v  that  we  have  calculated.  The  character- 

i-i.  c 

is tic  velocity  v£  was  calculated  using  the  relation 


v 

c 


(log  f  +  )vi4.VfJ_ 
Vi-L  +  Vf^ 


(19) 


11 


which  was  obtained  by  solving  eq.  (9)  for  v^.  The  characteristic  veloc¬ 
ity  determined  this  way  is  the  "experimentally  determined"  vc  because  the 
ion  fraction  f+  used  in  the  calculation  is  taken  from  experimental  data. 

It  can  be  seen  from  table  1  that  the  characteristic  velocity  v£  is  an  in¬ 
creasing  function  of  ip,  indicating  that  the  data  is  at  least  qualitative¬ 
ly  consistent  with  our  model  [see  eq.  (12)].  A  stronger  support  for  the 
model  is  indicated  in  fig.  2  where  log  v  is  plotted  against 

sin2<fr).  According  to  eq.  (12),  a  linear  relationship  is  ex¬ 
pected,  which  is  indeed  indicated  in  the  figure.  The  slope  of  the  line 
yields  the  parameter  -g-,  which  we  estimate  to  be  —  a  0.38.  We  note 
that  the  same  value  of  -g  was  obtained  with  the  1  keV  specular  data,  a 
strong  support  for  the  validity  of  the  model.  Assuming  that  the  poten¬ 
tial  energy  parameters  B  and  b  are  given  by  B  ■  3  keV  and  b  ■  5  &  *  as 

In  our  earlier  analysis  of  the  1  keV  data,  we  can  now  evaluate  zq  using 

A  A 

eq.  (11)  and  the  parameter  —  using  the  relation  —  *  v  exp(az  ). 

cL  &  C  O 

A 

The  calculated  values  of  z  and  —  are  also  shown  in  table  1.  The  param- 

o  a  r 

eter  —  represents  the  characteristic  velocity  at  z  *  0  and  is  charac¬ 
teristic  of  the  system  being  considered.  Its  values  calculated  at  differ¬ 
ent  'll  should  therefore  be  the  same  within  calculatlonal  and  experimental 

A 

uncertainties.  We  indeed  see  from  table  l  that  calculated  values  of  — 

a 

are  all  fairly  close  to  each  other,  lending  another  support  to  our  model. 

A  8  ^ 

Taking  the  average  we  obtain  —  ■  1.2  x  10  cm/sec  for  the  H  -graphite 

3l 

A 

system.  The  same  value  of  —  was  obtained  with  the  1  keV  data. 

With  the  parameters  —  ■  1.2  x  10®  cra/sec,  ■  0.38  and  B  »  3  keV, 

a  D 

we  now  can  calculate  the  characteristic  velocity  v£  using  eq.  (12).  The 
characteristic  velocity  calculated  using  eq.  (12)  represents  the  value 


predicted  by  our  model  and  thus  will  be  referred  to  as  "theoretically  de¬ 
termined"  vc»  In  table  2  and  fig.  3,  theoretically  determined  vc  are 
shown  along  with  experimentally  determined  v^  at  the  incident  energies  of 
1  keV  and  2  keV.  The  agreement  between  the  two  sets  of  values  is  seen  to 
be  very  good. 

We  now  look  at  the  nonspecular  data  of  Overbury  et  al.  In  table  3 

we  show  values  of  f+  and  v^  as  given  by  Overbury  et  al  along  with 

(or  9  -  <j»),  experimentally  determined  vc  and  theoretically  determined  v^ 

that  we  have  calculated  for  the  case  E^  (incident  ion  energy)  *  1  keV  and 

9  *  90°.  As  before  the  firjst  set  of  V£  (experimentally  determined  v^) 

was  calculated  using  eq.  (19)  with  the  values  of  f+  and  v^  reported  by 

Overbury  et  al.  The  second  set  of  v  (theoretically  determined  v  )  was 

A  & 

calculated  using  eq.  (15)  or  eq.  (18)  with  the  parameters  — ,  —  and  B 

<1  D 

determined  above  from  the  analysis  of  the  specular  data.  The  two  sets 

of  v£  given  in  table  3  are  shown  graphically  in  fig.  4.  Shown  also  in 

fig.  4  are  experimentally  and  theoretically  determined  v^  for  the  same 

scattering  angle  9  *  90°  but  for  a  different  incident  energy  *  2  keV. 

For  further  comparison  we  show  in  fig.  5  experimentally  and  theoretically 

determined  v  for  the  case  0  m  60°  and  E.  *  1  keV.  It  is  seen  that  the 
c  i 

characteristic  velocity  determined  by  our  model  is  generally  somewhat 
lower  than  the  experimentally  determined  v^.  The  general  qualitative  be¬ 
havior  of  the  two  sets  of  vc,  however,  coincides.  For  example,  both  ex¬ 
perimentally  determined  vc  and  theoretically  determined  v^  have  a  peak  at 
the  specular  configuration.  Considering  experimental  uncertainties  and 
uncertainties  in  scattering  cross  sections  and  system  parameters  (espec¬ 
ially  B  ■  3  keV  and  b  ■  5  X  \  which  are  rather  arbitrarily  chosen),  the 
agreement  is  good. 


I 

*S5 


3.2.  He  -Cu 

The  characteristic  velocity  of  the  He  -Cu  system  has  been  experimen¬ 
tally  Investigated  by  Bertrand  et  al  [12]  and  Verhey  et  al  [16,17], 

The  specular  data  of  Bertrand  et  al  [12]  shows  clearly  that  vc  Is 

0 

an  increasing  function  of  \f>  (*  -j)  at  a  given  ion  energy.  As  discussed  in 

the  previous  section,  our  model  also  predicts  the  increase  of  vc  with  i|/. 

/  i+u2  2  \a/b 

According  to  eq.  (12),  v  should  be  proportioned  to  (  — ^ —  sin  <jij  at 

a  given  energy  for  specular  reflection.  In  table  4  we  show  experimental 

values  of  v  that  we  have  estimated  from  the  data  of  Bertrand  et  al  as 
c 

well  as  theoretical  values  calculated  using  eq.  (12).  The  theoretical 

values  are  normalized  to  the  experimental  value  vc  =  2.0  x  107  cm/sec  at 

-  750  eV  and  ip  =*65  with  ^  ■  0.35,  the  value  which  seems  to  give  the 

best  fit  to  the  experimental  data.  One  can  observe  that  a  reasonable 

agreement  is  achieved  between  the  experimental  and  theoretical  values. 

The  value  of  —  can  now  be  estimated  by  requiring  that  it  yields  v  3 
cl  C 

2.0  x  107  cm/sec  at  E^  ■  750  eV  and  i|>  ■  65°,  If  we  take  B  *  3  keV,  b  =* 

5  X  *  following  Hagstrum  [8],  we  obtain  —  «  3.6  x  107  cm/sec.  If,  on  the 

3i 

other  hand,  we  choose  B  ■  1.805  keV  and  b  ■  3.867  8  *  according  to  the 

A  7 

table  given  by  Abrahamson  [21],  we  obtain  —  ■  3.0  x  10  cm/sec. 

2i 

Bertrand  et  al  have  also  provided  nonspecular  data,  which  indicate 
that  v^  increases  as  Che  incidence  angle  <|»  increases  from  grazing  angles 
toward  the  specular  condition  at  a  given  scattering  angle  9.  This  is  al¬ 
so  consistent  with  the  prediction  of  our  model.  Fig.  6  gives  the  plot  of 
the  characteristic  velocity  v£  calculated  using  eq.  (15)  for  different 
values  of  E^,  ip  and  9.  Two  different  choices  of  parameters  (^  ■  3.6  x 

107  cm/sec,  B  -  3  keV,  i  -  0.35;  and  -  -  3.0  x  107  cra/sec,  3  -  1.805  keV, 

b  a 

■  0.35)  yielded  virtually  the  same  values  of  v  ,  referred  Co  as 


theoretically  determined  V£  in  fig.  6.  Also  shown  are  experimentally 

determined  v  which  are  estimated  from  the  data  of  Bertrand  et  al.  Un- 
c 

fortunately,  only  relative  magnitudes  of  v£  can  be  calculated  from  their 
nonspecular  data.  We  therefore  have  normalized  experimentally  determined 
values  to  vc  =  1.36  x  10^  cra/sec  at  E^  *  750  eV,  i{t  ■  28.5°  and  0  *  90°. 
The  agreement  between  theory  and  experiment  is  not  very  good,  especially 
at  =  1250  eV.  This  is  mainly  due  to  the  fact  that  the  experimental 
data  does  not  show  any  visible  increase  in  v£  when  the  ion  energy  is  in¬ 
creased  from  750  eV  to  1250  eV,  where  such  behavior  is  inconsistent  not 
only  with  our  model  but  also  with  many  other  experimental  data.  We  note 
however  that  theory  and  experiment  do  agree  on  one  major  qualitative  fea¬ 
ture:  at  given  and  0  the  maximum  vc  occurs  at  the  specular  configura¬ 
tion.  We  also  note  that  both  theory  and  experiment  agree  that,  at  given 
and  <|>,  v£  increases  with  0  (this,  however,  is  not  clearly  indicated 
by  the  experimentally  determined  vc'at  *  1250  eV). 

Verhey  et  al  [17]  studied  the  He+-Cu  system  and  found  that  v is  an 
increasing  function  of  the  incident  ion  energy.  However,  it  is  difficult 
to  compare  their  analysis  with  ours  because  they  have  Included  collision- 
al  neutralization  and  ionization  of  the  atomic  type.  We  cau  still  obtain 

theoretical  values  of  v£  and  f+  using  eqs.  (12)  and  (9)  with  the  parame- 
A  ^  ^ 

ters  — ,  B  and  t-  for  He  -Cu  determined  from  the  data  of  Bertrand  et  al. 
a*  b 

Such  a  calculation,  however,  may  not  be  very  meaningful  because  these 
parameters  are  expected  to  vary  sensitively  with  surface  contamination. 
There  exists  only  a  remote  possibility  that  surfaces  used  by  two  differ¬ 
ent  groups  have  the  same  amount  of  surface  contamination. 

The  same  group  in  another  paper  [16]  has  reported  the  observation 
made  on  the  He+-Cu  system  that,  at  low  incident  energy  (E^  *  2  keV)  and 
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constant  0  (0  -  30°),  f+  is  approximately  a  symmetric  function  of  ij>  around 
the  specular  configuration  i|>  ■  j  0.  According  to  our  model,  a  complete 
symmetry  is  obtained  if  vf  “  v^»  if  u  ■  1.  For  He+-Cu,  the  func¬ 

tion  u  is  not  much  different  from  1  (u  =  0.992)  at  0  *  30°,  and  thus  an 
approximate  symmetry  is  expected,  as  was  observed. 

3.3.  He*-Ag,  Ne*-Ag  and  He+-Nl 

MacDonald  et  al  [19,  20,  22]  determined  the  characteristic  velocity 
for  the  systems  He+-Ni,  He+-Ag  and  Ne+-Ag  from  the  angular  distribution 
of  their  ion  yield  measured  at  constant  E^,  0  (=  90°)  and  <J»,  but  for  dif¬ 
ferent  exit  trajectories.  The  exit  trajectory  was  varied  by  measuring 
the  ion  yield  in  a  plane  perpendicular  to  the  plane  containing  the  inci¬ 
dent  beam  and  surface  normal  (the  incidence  plane).  The  ion  fraction 
measured  at  an  angle  $  to  the  incidence  plane  is  given  according  to  eq. 

(9)  by 


f 


+ 

* 


exp 


1 

v^sintji 


+ 


1 

VfCOS<j»COS$ 


(20) 


& 


If  the  variation  of  v^  with  $  is  neglected,  we  obtain 


i 

log  - - - £ —  f  — l—  -  i']  . 

f J,0o  V°8*  \  COS<£  ) 

(21) 


MacDonald  et  al  obtained  a  straight  line  when  log  — 1 —  is  plotted  against 

f*«o° 

v  cosi}i  (cosifr  ”0*  indicatlng  that  vc  varies  only  slightly  as  a  function 
of  $.  We  therefore  neglect  the  $  dependence  of  vc  in  our  analysis  to  follow. 
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MacDonald  and  O'Connor  [19]  obtained  the  energy  dependence  of  v  for 

c 


the  systems  He  -Ag  and  Ne  -Ag  by  evaluating  the  slope  of  the  log  — ^ — 

i  i  £**°° 

vs.  —  -1  line  at  different  energies.  In  table  5  we  show  Ef 

(scattered  ion  energy)  and  v  taken-  from  the  He+-Ag  data  of  MacDonald  and 

c 

A  *  * 

O'Connor  along  with  z  and  —  calculated  as  before.  Taking  the  average 

o  a 

A  7  + 

we  obtain  —  a  2.5  x  10  cm/sec  for  the  He  -Ag  system.  Fig.  7  gives  the 
plot  of  log  v^  as  a  function  of  log'  for  He+-Ag  and  also  for  Ne+-Ag. 

According  to  eq.  (12)  one  should  obtain  a  linear  relationship,  which  is 
approximately  indicated  in  the  figure.  The  slope  of  the  line  gives 

D 

From  fig.  7  we  estimate  ^  ■  0.35  for  He+-Ag  ahd  ^  ■  0.44  for  Ne+-Ag. 

3.4.  Ne~**-Au 

Brongersma  et  al  [15]  determined  the  ion  fraction  of  specularly  re¬ 
flected  particles  from  their  experimental  data  of  a  2  keV  Ne+-Au  colli¬ 
sion.  Shown  in  table  6  are  t|>  and  f+  taken  from  their  data  along  with 
experimentally  determined  characteristic  velocity  calculated  using  eq. 
(19).  We  again  clearly  see  that  v£  is  an  increasing  function  of  <»  -  9/2 
for  specular  configurations.  In  fig.  8  we  plot  log  v£  vs.  log  ( sin2* 
The  data  points  fall  approximately  on  a  straight  line  with  the  slope  esti- 
mated  to  be  ^  *  0.33.  Thus  we  see  that  our  model  is  supported  by 
the  experimental  data  even  for  the  case  of  heavy-ion  incidence.  Fr . lowing 

Hagstrum  [8]  we  choose  B  ■  36  keV  and  b  -  5  X  *  and  calculate  z  using 

o 

A  A 

eq.  (11)  and  —  using  the  relation  —  ■  vcexp(azQ).  The  calculated  val- 
A 

ues  of  z  and  —  are  also  shown  in  table  6.  The  values  of  z  calculated 

O  d  O 

are  seen  to  be  significantly  larger  than  the  values  0.03  -  0.3  X  quoted 


by  Brongersma  et  al,  possibly  dut  to  our  neglect  of  multiple  scattering 
which  may  be  significant  for  the  case  of  Ne+  incidence.  We  note,  however, 


that  the  values  of  —  calculated  at  different  <j»  remain  close  from  one 

& 


another,  with  the  average  -  *  5.85  x  10  cm/sec.  A  different  choice  of 

parameters,  B  3  13.625  keV  and  b  3  3.639  X  *,  as  suggested  by  Abrahamson 

[21],  was  found  to  yield  similar  values  of  z  and  somewhat  lower  values 

o 

A  A  7 

of  —  with  the  average  —  a  4.26  x  10  cm/sec.  The  two  different  sets  of 

dL  3 

the  parameters,  however,  yield  virtually  the  same  values  of  v£  which  are 
shown  in  fig.  9  along  with  experimentally  determined  v^.  The  agreement 
between  theory  and  experiment  is  surprisingly  good,  considering  that  our 
model  based  on  the  assumption  of  a  single  binary  collision  may  not  ac¬ 
curately  describe  the  Ne+-Au  collision  where  multiple  scattering  is  prob¬ 
ably  not  negligible. 

4.  Discussion 

The  characteristic  velocity  (more  specifically,  the  characteristic 
velocity  at  z  =»  0)  can  be  a  useful  parameter  for  surface  analysis  because 
it  is  sensitive  to  details  of  the  surface  such  as  the  molecular  composi¬ 
tion  and  the  electronic  state.  On  the  other  hand,  it  is  its  sensitivity 
to  surface  contamination  [13,23,24]  that  makes  both  experimental  work  and 
theoretical  analysis  difficult.  For  a  meaningful  determination  of  the 
characteristic  velocity,  the  surface  under  study  must  be  clean,  a  diffi¬ 
cult  experimental  task.  The  absolute  values  of  the  characteristic  veloc¬ 
ity  determined  by  our  analysis  may  not  represent  those  characteristic  of 
a  clean  surface.  MacDonald  and  Martin  [22]  reported  discrepancy  in  the 
value  of  the  characteristic  velocity  for  the  He+-Ni  system  determined  by 
two  different  methods  and  ascribed  it  to  the  contribution  from  colllsional 
neutralization  and  ionization  of  atomic  type.  There  however  seems  to 
exist  a  possibility  that  at  least  a  part  of  the  discrepancy  arises  from 
different  degrees  of  surface  contamination. 


The  main  purpose  of  the  present  work  is  to  test  the  validity  of  our 
model  introduced  in  ref.  18  and  further  developed  in  section  2  of  this 
paper.  It  is  encouraging  that  generally  good  to  excellent  agreements  are 
obtained  between  the  experimental  data  and  our  model,  as  described  in 
section  3.  This  Indicates  that  there  is  much  merit  in  our  model  accord¬ 
ing  to  which  surface  neutralization  is  a  dominant  mechanism  that  occurs 
between  low-energy  ions  and  a  solid  surface.  Atomic  neutralization  and 
ionization  and  a  host  of  other  processes  may  become  important  at  high 
energies,  but  in  the  low-energy  region  studied  here  (a  few  keV  or  lower), 
surface  neutralization  seems  to  be  the  key  process  that  determines  the 
ion  fraction  of  backscattered  particles. 
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are  calculated  as  described  in 


represent  experimentally  and  theoretically  determined 


0.025  0.0225  0.021  0.0215  0.0215  0.019  0.020  0.0185  0.0175  0.0195  0.018  0.020  0.019  0.0185  0.0195  0.018  0.0195  0.0175  0.016  0.0165  0.0155 


and  v  represent  experimentally  and  theoretically  determined 


E  =*  500  eV 


*  750  eV 


1250  eV 


^  =  10° 

ip  =  65° 

o 

o 

H 

n 

-9- 

i|>  ■  65° 

«p  -  10° 

»  65° 

v*(10^  cra/sec) 
c 

0.53 

1.7 

0.65 

2.0 

0.99 

2.0 

v^(10^  cm/sec) 

0.56 

1.7 

0.64 

2.0 

0.77 

2.4 

Table  4.  Characteristic  velocity  for  specular  configurations  for  the 
+  12 

He  -Cu  system,  and  v^  represent  experimentally  and 

2 

theoretically  determined  vc,  respectively.  vc  is  normalized 
to  the  experimental  value  of  2  x  10^  cm/sec  at  *  750  eV 
and  iji  a  65°. 
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ip(degrees) 

22.5 

30.5 

45 

53.4 

60.5 

69.1 

f+ 

0.009 

0.0132 

0.025 

0.025 

0.0266 

0.0322 

v  (107  cm/sec) 
c 

1.23 

1.48 

1.71 

1.91 

2.01 

2.02 

z  (X) 

0.97 

0.86 

0.74 

0.69 

0.66 

0.64 

o 

— (107  cm/sec) 

5.93 

6.00 

5.67 

5.88 

5.92 

5.70 

a 

Table  6.  Ne+-Au  data  for  specular  conditions  at  ■  2  keV. 

<J>  and  f+  values  are  taken  from  Brongersma  et  al  [15], 
and  v  ,  z  and  —  are  calculated  using  the  method 

CO  <1 

described  in  this  paper. 
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